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ABSTRACT 

We have carried out a survey of long 50ks XMM-Newton observations of a sample 
of bright, variable AGN. We found a distinctive energy dependence of the variability in 
RXJ0136. 9-3510 where the fractional variability increases from 0.3 to 2 keV, and then 
remains constant. This is in sharp contrast to other AGN where the X-ray variability 
is either flat or falling with energy, sometimes with a peak at ~ 2 keV superimposed 
on the overall trend. Intriguingly these unusual characteristics of the variability are 
shared by one other AGN, namely RE J1034-f 396, which is so far unique showing a 
significant X-ray QPO. In addition the broad band spectrum of RXJ0136. 9-3510 is 
also remarkably similar to that of RE J1034-I-396, being dominated by a huge soft 
excess in the EUV-soft X-ray bandpass. The bolometric luminosity of RX J0136.9- 
3510 gives an Eddington ratio of about 2.7 for a black hole mass (from the H beta 
line width) of 7.9 x IO^Mq. This mass is about a factor of 50 higher than that of RE 
J1034-I-396, making any QPO undetectable in this length of observation. Nonetheless, 
its X-ray spectral and variability similarities suggest that RE J 1034-1-396 is simply the 
closest representative of a new class of AGN spectra, representing the most extreme 
mass accretion rates. 
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1 INTRODUCTION 

Active galactic nuclei (AGN) are powered by gas accret- 
. ing onto the central super-massive black hole. As such their 
' intrinsic spectra are determined simply by the mass and 
spin of the black hole, together with the mass accretion 
rate, although their observational appearance is also af- 
fected by absorption. The discovery that the black hole mass 
scale s with the host galaxy bu lge properties l|Gebhardt et al.l 
I2OOOI : iMagorrian et al.l 1 19981 ) presented a means to clas- 
sify the bewildering variety of unabsorbed AGN spectra 
in terms of two of these three intrinsic properties. It was 
found that different optical spectral types clearly correlated 
with black hole mass and mass accretion rate, with the 
Narrow-line Seyfert 1 galaxies (NLSls: defined as Hp emis- 
sion line FWHMs le ss than 2000 km /s and [Oil I]/ Hp flux 
ratios less than 3.0 iGoodricbl 1 19891 ) as examples of lower 
mass black holes accreting at higher Eddington ratio than 
the typical broad line Seyfert Is (BLSls), which in turn 
had higher Eddington ratio/lower black hole masses than 
the LINERSs (e.g. B oroson 2002; Ostcrbro ck fc Pog gc 198^ ; 
Boiler. Brandt fc Finklil996l: lMathuj[200ol : iBian et 'aLil2008l : 



soft X-ray selected ROSA T AU Sky Survey l|Grupe et all 
1 19961 : iHasinger et alll 19971 ) . This strong ionising flux results 
in intense permitted optical Fe II emission lines from the 



to as coronal lines 



JWiUs et al.' 


1999: Kuraszkiewicz et al.l 


[2004: MuUanev & Ward 


20081). Other 



Winter et al.ll2009l ). 



NLSl's are intrinsically very luminous in the EUV/soft 
X-ray bandpass, comprising over half the AGN found in the 



observed chara c teristics include rapid X-r ay variability 
iLeighlvl Il999al : iBoller. Brandt fc FinkI 1 19961) and steeper 
soft X-ray spectra than thos e of BLSls ( Leighlvl Il999bl : 
iBrandt. Mathur fc Elvi3 1 19971 ). Taken together these dis- 
tinctive properties have led to them being consid- 
ered as the scaled up counterparts of the most lu- 
minous accretion states seen in stellar mass black 
hole binary (BHB) sys tems ie. those having " high" 
and "very high" states (IPounds. Done fc Osbornd Il995l : 
iMiddleton. Done fc Gierlinskill2007l v 

This link with the stellar mass black hole systems was 
strengthened by the discovery of an X -ray QPO in the NLSl 
RE J1034-I-396 (jGierliriski et al.ll2008l ). Such QPOs are com- 
mon ly strongest in the BHB observe d in their very high state 
(e.g. iRemUlard fc McClintockl l2006l ). Clearly it is important 
to understand whether QPOs are common amongst NLSl, 
or whether RE J1034+396 is unique in some other way. 

We first summarise the properties of RE J1034-)~396, in 
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order to provide the context for our study of RX J0136.9- 
3510, which is reported in this paper. The uniqueness of 
RE J1034+396 is suggested by its broad band spectral 
energy distribution (SED) which is extreme even amongst 
NLSl. It exhibits a peak in the far UV which connects 
smoothly onto the steep soft X-ray spectrum. These 
components form a huge "soft X-ray excess" with respect 
to the r ~ 2.2 X-ray t ail wh ic h dominates above ~ 2 keV 
"1995I: ICasebeer, Leighlv fc Baron 



JPuchnarewicz et al 



2006; M iddleton. Done fc Gierlinskil l2007i : iMiddleton et al 
2009i V The energy dependence of this variability is also 



very different to that commonly seen in other AGN. The 
fractional variability amplitude (as measured by root mean 
square, hereafter rms) rises steeply to ~ 2 keV and then 
levels off. This is most likely due to the presence of two sepa- 
rate components in the X-ray spectrum, with the variability 
being associated with the X-ray tail, whils t the soft excess 
comp onent remains more or less constant (|Middleton et al] 
I2OO9I ). This situation contrasts with the flat or falling 
rms spectra seen in other AGN, so metimes with a peak 
at 2 keV superimposed on this ( Vaughan et al.' '2003*; 
Vaughan et al. 20 04; Fabian et al . 2004j^ Gierlinski fc Done 



200e : iPonti et aLlbood : jPetrucci et ani2007l : iLarsson et all 
20081 ). which makes it more likely that the apparent soft 



excess in these objects is due instead to a single spectral 
component d istorted by reflection a nd/or absorption 
l|Crummv et al. 2006 ; Gierlinski fc Donell2006 ,') . 

In this paper we make use of the unusual energy depen- 
dence of the X-ray variability in RE J1034-I-396 to search for 
potentially similar objects. A survey of all long (> 50 ks) 
XMM-Newton observations of bright and variable AGN 
yielded a similar rms shape only in one object RX J0136.9- 
3510 (2MASSi J0136544-350952). This AGN also has a sim- 
ilar broadband spectrum to RE J1034-I-396, suggesting that 
they may form a subclass of the highest mass accretion rate 
AGN. 



2 SOURCE SELECTION 

We searched the XMM-Newton Master Log & Public Archive 
for pointed observations with exposure times in the PN in- 
strument of > 50000 s in "subject.category" "AGN, QSOs, 
BL-Lacs and XRB". This resulted in 115 observations. We 
further refined this criteria to include only those objects with 
PN count rates of > 1.0 counts/s, to include only bright 
sources for which the variability can be well determined. We 
also included NGC 4051 as this is a well known, bright, NLSl 
with an exposure of 49917 s (0157560101). Since this is not 
a statistical study its inclusion does not affect our general 
conclusions. 

We also searched the XMM-Newton Serendipitous 
Source Catalogue (2XMMi Version) for serendipious bright 
AGN detected in similarly long exposures, by setting 
pn_8Jlux ^10~^^ erg cm~^ s~^ and pn_ontime^ 50000 s. 
This yielded 29 observations. 

We then combined these two samples, and included only 
the longest exposure dataset or the best observation quality 
dataset for multiply observed objects. This gave a total of 
68 individual sources. 

For each source we calculated the fractional variability 
amplitude using Icstats in xronos 5.21. This is defined 



as the root mean square of the intrinsic (corrected for er- 
ror bars) variance about the lightcurve mean, normalised to 
this mean. Only sources which are strongly and significantly 
variable can provide constraints on the energy dependence of 
the variability. So we select only sources for which the rms is 
0.1 and where the uncertainty on this rms is ^ 0.33xrms, 
to ensure that the variability is significantly detected. This 
filtering leaves 15 sources. 

For each of these sources we calculated the fractional 
variaibilit y amplitude as a f unctio n of energy using the 
method of iGierhiiski fc Done! (120061 ). This showed the stan- 
dard range of rms spectral shapes. Only RX J0136. 9-3510 
(0303340101) displayed the very different type of rms spec- 
trum associated with the QPO AGN, RE J1034-(-396, in 
which the fractional amplitude rises as a function of energy, 
then remains high. Having identified this unusual AGN, we 
investigate its properties in more detail below. 



3 RX J0136-3510: LIGHTCURVE AND RMS 
SPECTRUM 

The lightcurve and rms spectrum are extracted from the X- 
ray data (0.3-10 keV), using SAS7.1.0 and xronos5.21. We 
use regions with radius of 16", 9.2" and 8.6" for pn, MOSl 
and M0S2, respectively. No significant pile-up is seen from 
the SAS command epatplot so no central region was ex- 
cluded. There are high background flares during the flrst 
10000 s, so we exclude these data, resulting in ~40000 s 
(with background count rate< 0.4 counts/s) for the rms 
spectrum generation. Figure [T] shows the resultant total 
lightcurve (PN, MOSl and M0S2), binned on 200 s. The 
source shows strong variability (also seen by Ghosh et al. 
2009, in preparation), with the rms fractional variation be- 
ing 0.13. However, there is no obvious QPO. 

We rebin the lightcurve on 2000 s in order to provide 
enough count statistics to calculate the fractional variability 
as a function of energy. Fig. [5] shows that this rms rises 
steeply from 0.3 keV to ~2 keV, then flattens off. The weak 
decrease from 2-10 keV is not signiflcant as the low count 
rate means their are large uncertainties on the last point. 
We also calculated the rms spectra from 200 s binning, but 
this is not significantly different. 



4 BLACK HOLE MASS ESTIMATION 

The opt ical spectrum has an Fell /H (3 fiux ratio of 
8.3 (Grupe et al. 19991'). compared to the ~ 1 average for 



NLSls (|Veron-Gettv. Veron fc Goncalvesll200]J), ma king RX 
J0146.9-3510 an unusual NLSl (|Ghosh et al.ll20o3 ).We use 



the optical spectrum shown in Fig [3] (D. Grupe, private 
communication) to estimate the black hole mas s. The Hfl 
line width is use d as a proxy for this (|Woo fc Urrvii2002 : 
iKaspi et al.ll2000l ) from 



Me 



4.817X 



- ALa(5100A) 0. 
lO-i^erys-i ' 



'^FWHM^ 



(1) 



Comparing this method for a sample of AGN with re- 
verberation mappi ng, the rms difference is about 0.5 dex 
(|Woo fc Urrvll2002l '). The flux at A = 5100A from the optical 
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Figure 1. Background subtracted lightcurve of RX J0136. 9-3510 
binned on 200 s. The exposure start time (UTC) is 2005-12-14 
20:45:30, but the first 10 ks was excluded due to the high back- 
ground contamination. 
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Figure 2. The rms spectrum correspongding to Figure [T] the 
binning time is 2000s. 
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Figure 3. The optical specturm of RX J0136. 9-3510, including 
the H/3 emission line fitting. The data is from Grupe D. and we 
show this figure with his permission. 



Using our new value for the FWHM, the resultant black hole 
mass is 

Even with this lower black hole mass, this AGN is still 
probably ~ 10 x more massive than the QPO AGN, RE 
J1034-I-396, (with the previous mass estimate being ~ 50x 
larger). Thus any similar QPO in RX J0136. 9-3510 would be 
on timescales 10-50 x larger, requiring a much longer X-ray 
observation in order to detect it. 



spectrum is 
tance is Dl 



-1.5x10 
1455.2 Mpc for z 



1 fi — 9 - 

erg cm s 



The luminosity dis- 
0.289 assuming Ho = 72 
km Mpc~\ Q,M = 0.27, and Q,vac = 0.73. 

Estimating the FWHM of the line is not straightfor- 
ward as the H/3 line profile is complex. There is clearly also 
a component from the extended narrow line region (NLR). 
This should be similar to the profile of the [OHI] emis- 
sion line. So lGrup^ et al.l (|l999l ) use a template constructed 
from the [OIIIJA5007 line to represent the narrow compo- 
nent, together with a broader Gaussian component, with 
FWHM=1320 km/s, to reconstruct the H/3 line. This gives 
a black hole mass estimate of 1.3x lO'^ Mq. 

However, the optical spectrum plainly has limited 
signal-to-noise, and the [OHI] 5007 line is probably contam- 
inated by Fell emission as it should be in 3:1 ratio with the 
(unobserved) [OHI] line at 4959A. Hence we perform our 
own best fit of the H/3 fine profile with a gaussian of width 
870 km/s for the narrow line component, see Figure O This 
gives a FWHM for the broad component of 3200±2600 km 
s~^, with a (not significant) blueshift of -370±1100 km s~^. 



5 BROAD BAND SED ANALYSIS AND 
EDDINGTON RATIO 

We use the standard products to obtain the XMM-Newton 
X-ray (PN) and optical/UV (OM) spectra. We then combine 
these with selected continuum points from the (non simulta- 
neous) optical spectum (D. Grupe, private communication) 
using FLX2XSP to incorporate these into the same format as 
the XMM-Newton data. We likewise include the J,H and 
K near infrared fiux points from 2MASS, and perform all 

spectral fitting using xspecll . 3.2. 

We follow the approach of IVasudevan fc FabianI (|2009l ) 
in modelling the broadband SED, using diskpn to model 
an accretion disc extending down to the last stable orbit 
around a non-spinning black hole. However, our source is 
at 2 = 0.289 so we modify the diskpn code to incorpo- 
rate the redshift dependance. And also we use the higher 
estimated black hole mass in Section |3] to calculate the 
normalization of the diskpn and fix it at 2910 during the 
spectral fitting. Compton scattering of these disc photons 
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Table 1. The best-fit parameters in our model for the optical to X-ray broad band SED. 
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' •'''^-0.31 


n 00+0.03 ir, 17+0.72 
^•^^-0.02 ^^--^'-0.75 




1 10+0.55 
^•^^-0.37 


475/493 



can be approximated by a broken power law (bknpower), 
with index of F = 0.33 below a break at SkTdisk- We 
then added a low temperature Comptonisation component 
to model the soft excess, using comptt, with seed pho- 
tons set to the disc temperature. We assume that these in- 
trinsic components are absorbed by both gas and dust in 
our Galaxy, and so fix this parameter to the value of the 
Galactic gas columi£] (wabs) of Nh = 0.0208x10^^ cm'^. 
The reddening (red den) is linked to this a ssuming a E(B- 
V)=1.74xlO"^^ArH (|Spitzer fc Lvmanll978l ). Since the opti- 
cal data were not simultaneous with the UV and X-ray data, 
we allow for long time-scale variation as a constant offset in 
normalisation between the XMM-Newton data and the op- 
tical spectrum. We exclude the infrared data points from 
our spectral fitting, since the model we use describes the in- 
trinsic emission from the accretion flow whereas the infrared 
emission is likely due to reprocessing of the UV emission by 
dust in the host galaxy, plus a possible contribution from in- 
trinsic starlight. The resultant best-fit parameters are given 
in Table [U 

Figure |3] shows the rebinned data (black), with the 
optical flux corrected for their best flt normalisation of 
~ 0.92 X that of the XMM-Newton UV and X-ray spectra, 
and all datapoints are corrected for absorption/reddening. 
This model is a good description of the overall shape of the 
optical/UV/X-ray spectrum and gives a bolometric (0.001- 
100 keV) flux of l.lx 10~^''erg cm~^s~^, corresponding to a 
bolometric luminosity of 2. 7x10*^ erg s~^. 

Using our own fltting value for the FWHM, from which 
the resultant black hole mass is 7.85x10'' Mq, the Edding- 
ton ratio for RX J0136. 9-3510 is as high as 2.7! Moreover, 
if we adopt the lower estimated black hole mass which is 
1.3x IO'Mq, then the normalization of diskpn is 80, and the 
resultant disk temperature rises to 31 eV. The UV region is 
then dominated by disk emission, though the Comptonisa- 
tion still is required to model the soft X-ray excess. This 
gives a bolometric flux of 8.8 x lO"'^^ ergs s~^ cm~^, giving 
an even higher Eddington ratio of 13. 2! This is the highest 
known Eddington ratio for an AGN l|Vasudevan fc FabianI 
I2OO9I : IShen et all lioosl ). These two best flt models basing 
on different black hole masses set a range for the Edding- 
ton ratio of RX J0136. 9-3510, making it a AGN with super 
Eddington accretion flow. 



6 SUMMARY AND CONCLUSION 

In this paper we report that RX J0136. 9-3510 is the only 
well observed. X-ray bright, variable AGN which has a sim- 
ilar energy dependence to its X-ray variability as the so far 
unique QPO AGN, RE J1034-I-396. Its Eddington ratio is 
similarly high, around 3, although its larger mass means 
any QPO is undetectable in our data. Its broad band SED 

^ http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 
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Figure 4. RX J0136. 9-3510 unfolded spectrum. All data points 
from diff'crent wave bands are included in this figure, though the 
infrared points are not included in the model fitting. The model 
spectrum without galactic extinction and dust reddening is also 
generated and superposed on the source spectrum, with red rep- 
resenting redshifted disk component, green representing compTT 
and blue representing bknpower. The orange line shows the total 
model spectrum. 



is also remarkably similar to that of RE J1034-I-396, being 
well modelled by a low temperature, optically thick Comp- 
tonisation of the accretion disc spectrum, plus a tail extend- 
ing to higher energies. S pectra such as this have also been 
fit by "slim" disc models (|Abramowicz. Kato fc Matsumotol 
1 19891 ) ■ where the accretion rate is so high that radiation 
cannot easily escape vertically b efore it is carried radially 
(advected) along with the flow (|Puchnarewicz et al.ll200"ll : 
I Wang fc Netzeij |2003| ). However, simple slim disc models 
do not fit the curv ature of the soft X-ray spectra as well 
as Comptonisation l|Middleton et al]|2009l ). although more 
complex models of slim discs do include s uch scattering in 
the disc atmosphere (e.g. lKawaguchill2003l ). 

Low temperature, optically thick Comptonisation is 
also occasionally seen in the stellar mass black hole bi- 
nary systems, for example in the most e xtreme mass ac- 
cretion rate spectra of G RS 19154-105 (jMiddleton et al] 
l2006l : lMiddleton et al.|[2009l ) . Recent studies of spectra of the 
Ultra-Luminous X-ray sources also indicates that these are 
well m odelled by such material (|Gladstone. Roberts fc Done! 
I2OO9I ). This evidence suggests that there is indeed a dis- 
tinct spectral st ate which can only be attained b y super 
Eddington fiows l|Gladstone. Roberts fc Don3l2009l ) . Future 
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long duration X-ray observations of AGN should reveal ad- 
ditional examples, and in those cases with low black hole 
masses, the presence of a QPO is a prediction. 
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